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Abstract 

Inclusive differential cross sections do p A/dxF and du p A/dpl for the production of Kg, A, and A particles 
are measured at HERA in proton-induced reactions on C, Al, Ti, and W targets. The incident beam energy 
is 920 GcV, corresponding to ^fs = 41.6 GeV in the proton-nucleon system. The ratios of differential cross 
sections d<jpA(K-%) / do p A(A) and do p A{A) / du p A{A) are measured to be 6.2±0.5 and 0.66±0.07, respectively, 
for xpK— 0.06. No significant dependence upon the target material is observed. Within errors, the slopes 
of the transverse momentum distributions du p Ajdp\ also show no significant dependence upon the target 
material. The dependence of the extrapolated total cross sections <j p a on the atomic mass A of the target 
material is discussed, and the deduced cross sections per nucleon a pN are compared with results obtained at 
other energies. 



Key words: Hypcrons, cross section, A-dcpendence 
PACS: 13.85.-t 13.85.Ni 13.85.Lg 



1 Introduction 



We present measurements of inclusive production cross sections for K s , A, and A particles, collectively 
referred to as V° particles, in collisions of 920 GeV protons with several nuclear targets of different atomic 
mass A. The production of V° particles has been studied by numerous experiments with different beams Ok, 
K, p, n, £~) on different targets, covering a momentum range of 4-400 GeV/ c (see (|, |[ ||, [| ||, |[ fj], ||, [| 
|lH O, [l3| [l4|, [ll| E^, [l7| [IJ, [H], ^(J and refs. therein). A large fraction of these experiments utilized bubble 
chambers. At center-of-mass energies above yfs « 30 GeV, inclusive V° production cross sections have been 
measured by experiments at the CERN proton-proton Intersecting Storage Rings (ISR) |2^, 
However, only a lower limit for the A cross section has been reported ^2|. Also, the V° total production 
cross sections reported in [ p4[ are substantially above expectations based on extrapolation of results obtained 
at lower energies. 

Recently, progress in heavy ion physics has renewed interest in studies of strange-particle production. 
One of the main goals of heavy-ion experiments is the observation of the quark-gluon plasma 1 26 , and one 
of the signatures for this new state is enhanced production of strange particles J2?J . Observables of interest 
are ratios of antibaryons to baryons at mid-rapidity, which are important for net baryon density evaluations 
and which have been used recently to extract values of chemical potentials and temperatures |2^, |2^] for 
Au-Au collisions at RHIC. Also of interest are the squared transverse momentum {pi) distributions, da /dpi, 
which contain information about the temperature of the system; specifically, if it reached thermal equilibrium 
(see, e.g., pp|). These investigations attach importance to a comprehensive measurement of strange-particle 
production properties in "ordinary" nucleon-nucleon (NN) and nucleon- nucleus (NA) collisions. The latter 
are expected to establish a valuable baseline for comparisons among AA results p8] . 



2 The HERA-B Experiment 

HERA-B is a fixed target experiment at the 920 GeV proton storage ring of HERA at DESY It was 
originally designed to study CP violation in B° — ► J/ipKg decays; thus, the spectrometer was optimized 
to detect and precisely locate decays of J/ip and longer-lived Kg mesons. The first experience with the 
apparatus was gained during an engineering run in the year 2000. While the main goal was to commission 
the sophisticated J/^-based trigger system, there were also run periods dedicated to recording minimum 
bias events, i.e., events selected with a random trigger that uniformly sampled all HERA bunches. These 
data allow for precise measurements of strange-particle production cross sections. 

A plan view of the HERA-B spectrometer is shown in Fig. [j]. The spectrometer dipole magnet provides 
a field integral of 2.13 T-m, with the main component perpendicular to the x-z plane. The apparatus 
(including particle identification counters) has a forward acceptance of 10-220 mrad in the bending plane 
and 10-160 mrad in the non-bending plane. The experiment uses a multi-wire fixed target which operates 
in the halo of the proton beam during HERA e-p collider operation. Up to eight different targets can be 
operated simultaneously, with their positions being adjusted dynamically in order to maintain a constant 
interaction rate between 1 and 40 MHz. For the measurements presented here, wires made of carbon (C), 
aluminum (Al), titanium (Ti), and tungsten (W) were used. Their transverse width was 50 /an, and their 
thickness along the beam direction was 1000 /im for carbon and 500 /im otherwise. 

The tracking system consists of a vertex detector system (VDS) and a main tracker system. The latter is 
separated into an inner tracker (ITR) close to the proton beam-pipe and an outer tracker (OTR) farther out. 
The VDS |32) features 64 double-sided silicon microstrip detectors arranged in eight stations along, and four 
quadrants around, the proton beam. The silicon strips have a readout pitch of approximately 50 fiva.. The 
strips are rotated on the wafer such that each pair of wafers provides four stereo views: ±2.5° and 90° ±2.5°. 
The ITR uses GEM microstrip gas chambers |}3) , and the OTR uses honeycomb drift chambers |54| . There 
are a total of 13 ITR + OTR tracking stations, with each one referred to as a "superlayer." For the analysis 
presented here, only data from the six superlayers located downstream of the magnet are used. Also, this 
analysis does not use particle-identification information from the ring- imaging Cherenkov counter Pq| , the 
electromagnetic calorimeter j36|, or the muon detector ]37[| . 



3 Data Analysis 



The results presented here are based on a sample of ~ 2.4 million randomly triggered events recorded in 
April 2000. The HERA proton beam had 180 filled bunches and a bunch crossing rate of 10.4 MHz. The 
interaction rate was adjusted to values between 2 and 6 MHz, leading to most bunch crossings (> 90%) 
having either or 1 interaction. For this analysis, only tracks with a minimum of five hits in the VDS and 
ten hits in the main tracker were used. The same event selection procedures were applied as were used to 
determine the integrated luminosity. No explicit selection on track multiplicity was applied. 

In each event with at least two tracks, a full combinatorial search for V° candidates was performed. V° 
candidates were selected from all pairs of oppositely charged tracks that formed a vertex downstream of the 
primary vertex. The minimum distance between the two tracks was required to be less than 0.7 mm. No 
particle identification criteria were applied to the tracks. The primary vertex was determined from all tracks 
reconstructed in the VDS excluding the V° tracks. Using tracks reconstructed in both the VDS and the 
main tracker, we measured the spatial resolution of the reconstructed primary vertex to be 0.7 mm along 
the z-direction. The position of the primary vertex was required to coincide with the center of a target 
wire to within three standard deviations (2.1 mm). If a primary vertex could not be reconstructed, the z 
coordinate of the target was used to calculate the ^-component of the particle's flight length. An Armenteros- 
Podolanski plot j3j| for the selected V° candidates is shown in Fig. |2|. Clusters of events shaped according 
to the kinematics of K s -^ir + Tr~, A^pir~ and A^pTT + decays are clearly visible. The V° candidates 
from those regions of Fig. in which two of the three V° species overlap -and thus are kinematically 
indistinguishable - were rejected. This reduced the K s yield by 3.5% and the A (A) yield by 10% (20%). 
An additional requirement pt • ct > 0.05 (GeV/c) • cm was applied to reduce background from 7—* e + e _ 
conversions, where p t is the transverse momentum relative to the V° line-of-flight and r is the V° proper 
lifetime. This requirement also reduced combinatorial background from the target, which populates the 
lower region of Fig. |[ 

The invariant mass distributions for the selected candidates are shown in Fig. ||. Clear signals correspond- 
ing to Kg, A, and A particles are visible. In order to estimate the event yields, fits with a Gaussian function 
for signal and a third-order Legendre polynomial for background were performed. The width obtained for the 
Kg mass peak was about 5.6 MeV/c 2 , which is consistent with the momentum resolution oia(p)/p 2 10~ 4 
(GeV/c)" 1 . The resultant event yields are summarized in Table [j] for each target material. 



4 Acceptance Determination 

A Monte Carlo (MC) simulation is used to determine the reconstruction efficiency for the selected particles 
and decay channels. The production of a V° particle in an inelastic event is simulated using a kinematic 
distribution of the form: 

^-=C(l-\ XF \ r .eM-B^ t ), (1) 

where if is the Feynman-ai variable and pt the transverse momentum in the laboratory system. This 
phcnomenological ansatz is motivated by quark counting rules and phase space arguments [|40| and has been 
shown to describe a substantial amount of data well (see [^0| and refs. therein). The value of the parameter 
B of 2.1 (GeV/c) -2 is taken from Ref. {MJ, and a flat xf distribution is assumed (n — 0). The xf bin size of 
0.015 is chosen to be compatible with the momentum resolution. After the generation of the V° particle, the 
remaining momentum is assigned to a virtual 7r + , which is input as a beam particle to the FRITIOF 7.02 
package [Q in order to further simulate interactions within the nucleus. 

The generated particles are propagated through the geometry and material description of the detector 
using the GEANT 3.21 package Q. The detector response is simulated including realistic descriptions of 
chamber efficiencies and dead channels. The MC events are subjected to the same reconstruction chain as 
that used for the data. For tracks within the geometrical acceptance that originate from V° decays, we 
determine an efficiency of ~ 90% to reconstruct the track and assign a momentum. This value is confirmed 
by the analysis of real data. 

After all cuts, the efficiency to identify a V° particle is approximately 10% for K s and 5% for A and A. 
These values are dominated by the geometrical acceptance. All efficiencies are computed in bins of rapidity 
y and pf in order to be independent of the details of the MC model of V° production. 



5 Production Cross Sections 



The cross section Aayo for the production of a V° particle within the spectrometer acceptance can be 
expressed as: 

1 ff N v „(y,p 2 t ) 



Aa v° = JT^ / / ( 2 d V ■ d Pt . (2) 

b v « £ JJaccp £\y>Pt) 

where N v o(y,p 2 ) is the number of V° candidates observed in bins of rapidity y and p\. The detection 
efficiency e(y,Pt) is calculated from the MC simulations described above. The branching fraction Bya for 
the detected decay V° — > / is taken from Ref . |4j| . Finally, the integrated luminosity £ is computed from the 
interaction rate at the target |38| . The resultant luminosities for each target material are listed in Tabic [j]. 

The efficiency-corrected values of the inclusive differential cross section d<j p A/d,XF for V° production 
within the spectrometer acceptance are listed in Table |]for the various target materials. The corresponding 
xp interval is [—0.12,0], except in the case of Kg production on the W target, where it is [—0.09,0]. 

As the x f acceptance of the spectrometer is restricted, for comparison with results from other experiments 
it is necessary to extrapolate the differential cross sections da p A/dxF to the entire kinematic range xp G 
[— 1, +1 ]. This was done using the parameterization da/dxp oc (1— |x^|) n , where n is a constant. The values 
used for n are taken from previous measurements of inclusive strange-particle production p0| and are listed 
in Table ^. The difference between values for different V° species is compatible with theoretical predictions 
based upon quark counting rules Q. The resulting V° total production cross sections {<t p a) are also listed 
in Table |^. The fraction of the cross sections within the detector acceptance is 30% for K s (n = 6.0), 17% 
for A (n = 2.2), and 34% for A (n = 8.0). 

The dependence of the measured cross sections on the atomic mass of the target material is shown in 
Fig. ^ along with fits of two different kinds. The dashed lines are fits to the form o p a oc A a , while the 
solid lines are fits within the framework of the Glauber model pl| . These latter fits allow us to extract 
the production cross section per nucleon, <j p n- All fit results are listed in Table |[ The Glauber model 
calculations use a nucleon density given by the Saxon- Woods potential ftiql , and the total cross sections for 
KN and AiV collisions and the slope of the diffraction scattering cone are taken from Refs. (4^, [43]] . 

5.1 The pf Differential Cross Sections 

The extrapolation of the the pi distributions for Kg, A, and A particles over the entire xf range [— 1, +1 
yields the differential cross sections do jdp\ listed in Table |[ These data are presented graphically in Fig. 
together with fits to the form 

— I = <T-B-exp(-B-pt), (3) 
a Pt 

where the parameter B is independent of 

andp?. Table | summarizes the values of B obtained for the 
different target materials and V° particles. No significant dependence upon the target material is seen. 



5.2 Systematic Uncertainties 

To estimate the uncertainty in the efficiency determination, we varied the V selection criteria applied to data 
and simulated events. Different track quality requirements used in the vertex reconstruction were studied in 
order to understand differences observed between data and MC results; the variation of cuts causes a change 
of 7% in the cross section values. Cuts on the primary and V° vertices give a 2% change. The variation 
of the distance required between tracks to form a V° decay vertex contributes another 2%. The systematic 
error due to the limited statistics of MC events is 1% for Kg and 3% for A and A particles. The use of 
different fitting functions for the shape of the background in the invariant mass spectra, and the variations 
of bin sizes for invariant mass, xf, and pi, result in a change in the cross sections of approximately 5%. 

The efficiency of the ITR improved during data taking (as this system was further commissioned), and 
this caused a systematic shift in the track reconstruction efficiencies. These shifts are +14%, +11%, and 
+20% for K s , A, and A particles, respectively. 

We used three independent methods to determine the integrated luminosity. The results differ by 4-7%, 
which we include as an additional systematic uncertainty. For more details, see Ref. [p8[. 



To extrapolate do /dpi to the entire kinematic range of Xf , we used Eq. |l] with experimentally determined 
values pOfl of the parameter n. Varying this parameter by the estimated experimental uncertainty (±0.5) 
leads to a variation of 5% in the Kg cross sections, 12% in the A cross sections, and 2% in the A cross 
sections. 

The cross sections were measured on different nuclei, and extracting the cross sections per nucleon (o p n) 
from the measured values was done via the Glauber model. Variations of the parameters of this model give 
a 12% systematic uncertainty in <j p n- 

Adding all systematic errors in quadrature, we obtain the total systematic errors for the cross sections 
per nucleus do p A/ 'dxp, da p A/dp\, and <j p a, and for the cross section per nucleon o p n • All systematic errors 
are given in Tables ||, and ||| (listed after statistical errors) ; the errors are asymmetric due to changes in 
running conditions for the ITR. 



6 Discussion of Results 

More precise than the production cross sections themselves are cross section ratios, since, for these, acceptance 
corrections and systematic errors to a large extent cancel. From the measurements presented above we 
calculate the ratios da{K s )/da{A) = 6.2 ± 0.5 and da (A) /da (A) = 0.66 ± 0.07, for x F « -0.06 (i.e., 
mid-rapidity). The da(A)/da(A) ratio is plotted in Fig.^J for the different target materials; no significant 
dependence upon the target material (i.e., atomic mass A) is seen, and this holds also for the da(K s )/da(A) 
ratio. 

Ratios between antiparticle yields and particle yields measured at mid-rapidity have attracted special 
interest in studies of the dependence of AA collision dynamics upon collision energy p8j, p9| . Fig.|6| shows 
the ratio of yields for A and A particles at mid-rapidity in AA collisions as a function of center-of-mass 
energy ^/snn- The data plotted are from the BNL AGS (JsWn ~ 5.4 GeV), the CERN SPS (^saw- 10- 
30 GeV), and RHIC (y^Tjv = 130 GeV ) p[ ||, pj |g, |50] |§. Comparing these measurements 



with corresponding results from pp [TtJ |46|| and pA |}46|, \H[ 48j] collisions shows these data to exhibit a 
similar energy dependence, and the HERA-B result at a/s = 41.6 GeV is consistent with this trend. Closer 
inspections suggests that the pp, pA, and AA data follow individual curves that are slightly shifted in energy. 
This would be consistent with results from NA49 which found that in pp collisions at y/s nn = 17.3 GeV 
the ratio da (A) /da (A) is about 60% larger than that in pA collisions, and in A A collisions at the same ^/snn 
the ratio is about 60% smaller. However, more data are needed to confirm this behavior. 

The measured p\ spectra presented above are well-parametrized by Eq. |[ The results (Table ||) show no 
significant dependence upon the target material, and they are compatible with results from other experiments, 
e.g. from pp collisions at ^/s = 27.6 GeV, where at pj < 1.5 (GeV/c) 2 B parameters for Kg, A and A particles 
of 3.59±0.18 (GeV/cT 2 , 2.73±0.21 (GeV/c)~ 2 and 2.66±0.74 (GeV/cT 2 , respectively, have been extracted 
p8[ ; for more examples see the compilation in j^O). Due to the limited statistics available, the present study 
was done in the p\ region below 1.1 (GeV/c) 2 ; thus we cannot confirm the flattening of the p\ spectrum 
above 1.2 (GeV/c) 2 (see e.g. @). 

The narrow acceptance in xf for this data set also precludes us from measuring "leading particle" effects, 
which have been observed by other fixed target experiments at xf > 0.4 (see e.g. |2(], [53|]). However, we 
have measured for the first time V° differential cross sections for different target materials at negative xf 
(—0.12 < xf < 0). We extrapolate these results to the full xf range to obtain total cross sections for different 
values of atomic mass A. We fit these values to the conventional expression a P A oc A a and also to the 
Glauber model. The values obtained for a from the first fit (Table |3|) are very similar to those observed in 
other hadroproduction processes, which can be characterized J|4| by a(xp) — 0.8 — 0.75xf + 0.45 x f /\xf\ ■ 

From the Glauber- model fits we extract V° cross sections per nucleon, a V N (Table |3|). Fig. |7| presents 
a comparison of these results with previous experimental results [|[ 0, ||, S, [l(], |ll|, O, [U|, |IJ [l5|, [l6| |T^, 
[jj, II, [p, HH ||. The HERA-B results show good agreement with the general systematic trend as 



quantified by the curves. These curves were calculated from a fit [Q for the average number of V° particles 
per inelastic collision as a function of ln(s), in the interval 13.5<i/s<28 GeV. The data in Fig. [7 suggest 
that this parametrization is valid up to the HERA-B energy of y/s = 41.6 GeV, and even up to the highest 
ISR energy of ^fs — 63 GeV. Our value for the Kg cross section is in good agreement with the value obtained 
at v / s = 52.5 GeV p3[|. The large cross section reported in p4[ for y's = 63 GeV has been recognized |2j| 



to be biased by ghost tracks. In addition, our value for the A total cross section per nucleon is in good 
agreement with the values obtained at the ISR at v / s = 53 GeV and y / s = 62 GeV [^2|. As A production is 
known to receive large contributions from fragmentation processes, such a good agreement is noteworthy: it 
implies that also for small xf the parameterization of the (double-differential) cross section given in Eq. |l| 
is valid. 
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Table 1: Integrated luminosities Ha in mb -1 for the indicated targets of atomic mass A [[58], and the 
corresponding numbers of events N(evt) and V° particles N(V°) obtained. Uncertainties listed for La are 
systematic; those listed for N(V°) are statistical. 





C 


Al 


Ti 


W 


A 


12 


27 


48 


184 


£A-mb 


1093 ± 38 


1030 ± 67 


308 ± 10 


47 ± 3 


N(evt) 


496694 


893885 


467943 


512897 


N(K° S ) 


2566 ± 58 


4398 ± 81 


2022 ± 52 


909 ± 40 


N(A) 


512 ±31 


831 ± 43 


412 ±27 


224 ± 24 


N(A) 


241 ± 24 


503 ± 41 


222 ± 23 


144 ± 22 



Table 2: The inclusive differential cross sections da p A/dxF in mb measured in the Xf interval [—0.12, 0] for 
the production of Kg, A, and A particles on indicated targets. Also listed are the total cross sections <j v a 
obtained by extrapolation to the full xf range (see text). Values are followed by statistical and systematic 
errors. The values of the parameter n used for the extrapolations (see Eq. [I]) are also given. The xf interval 
for the case of K s production on the W target is [—0.09, 0]. 





Trgt 


n 


da p A/dxF (nib) 


(T pA (mb) 




C 
Al 


6.0 


215. ± 12 + f 3 
429. ± 33 + If 


86. ± 5 + f 
174. ± 13 + l° 6 




Ti 




612. ± 33 + 


248. ± 13+ 




W 




2044. ± 167.±fl|; 


761. ±63 + 1 !!; 


A 


C 


2.2 


53. ± 4+ ' 3 \ 


38. ± 3.1 




Al 




90. ± 8+ 1 !; 


64. ± 6.1^. 




Ti 




163. ± 13.+^; 


115. ± 20+^ 




W 




558. ±62 .tit 


399. ± 47.1^ 


A 


C 


8.0 


32. ± 3+ ' 2 


11. ± 1.1 ™ 




Al 




69. ± 8+ 1 !; 


25. ± 3 + |;f 




Ti 




96. ± ll.± 2 6.- 


34. ± 4 + i 5 7 




W 




375. ±58.±||; 


131. ±20+^; 



Table 3: The V° production cross sections per nucleon (ct p at) in mb with statistical and systematic errors, 
and the values of a resulting from fitting the cross sections per nucleus to the form <j v a oc A a . For the values 
of a, the errors listed are statistical only. 





K° 
A S 


A 


A 


a pN (mb) 


9.0 ± 0-3±i:| 


4.3 ±0.2!°;? 


1.4 ± 0.1±g:| 


a 


0.78 ±0.04 


0.85 ±0.05 


0.82 ±0.07 



Table 4: The inclusive differential cross sections da v Ajdp\ in mb/ (GeV/ c) 2 for the production of V° particles 
on indicated targets as obtained from extrapolation to the full xf range. Values are followed by statistical 
and systematic errors. The p 2 bins (Ap 2 ) are in (GeV/ c) 2 ; the p t resolution is 1.7 MeV/c. Because of limited 
statistics, no cross sections are given for the production of A particles on the W target. 



Ap 2 


C 


Al 


Ti 


W 




0.0- 


0.2 


230. ± 20 +H, 


454. ± 54 +11; 


619. ± 53+^; 


1901. ± 254 +^f; 


0.2- 


0.4 


94. ± lO.t 1 !; 


191. ± 28 +f 7 


275. ± 29.1^; 


899. ± 141.1 1 |f' 


0.4- 


0.6 


49. ± 7.1 I] 


101. ± 19 + 1 g; 


156. ± 20+12; 


560. ± 92 + H\ 


0.6- 


0.8 


28. ± 5.1 i 


64. ± 14 + 1 *; 


99. ± 15.1^; 


229. ± 57 + 21; 


0.8- 


1.0 


18. ± 4 + f 


33. ± 10.1 3; 


51. ±11 + 1; 


146. ± 47.1 2 |; 


1.0- 


1.2 


14. ± 3 + 2 ; 


29. ± 10.1 3. 


42. ± 10.1 3; 


75. ± 34 + ^; 


A 


0.0- 


0.2 


58. ± ll.t 1 "; 


103. ± 23.1 2 ^; 


91. ±34+^; 


677. ± 200 + 1 g^; 


0.2- 


0.4 


53. ± 6.1 ?; 


78. ± ii-±ii; 


167. ± 20+f 2 : 


327. ± 81.1 %■ 


0.4- 


0.6 


32. ± 4 + |; 


54. ± S.l^ 


94. ± 13. li| 


364. ± 60.1 H: 


0.6- 


0.8 


17. ± 3 + | 


33. ± 5 + |; 


70. ± 10. +_ l l\ 


249. ± 45.1 35; 


0.8- 


1.0 


15. ± 3.1 l[ 


27. ± 5.1 l\ 


57. ± 10. 1 1 "; 


263. ± 45.1 f 7 : 


1.0- 


1.2 


8.± 2+ i. 


14. ± 3.1 l\ 


44. ± 8.1 l 


120. ± 31.1 2 f; 


A 


0.0- 


0.2 


19. ± 5. 1?; 


49. ± 11. 1"; 


71. ± 17.1 1 ! 




0.2- 


0.4 


16. ± 2. 1?; 


29. ± 4.1 I] 


29. ± 6.1 2. 




0.4- 


0.6 


8- ± i. ±S:"g 


16. ± 3.1 1; 


28. ± 4 + 6 2 




0.6- 


0.8 


s- ± i. ±S: 3 


i6. ± 2.1 1: 


20. ± 4.1 2. 




0.8- 


1.0 


4.6 ±0.9tS; 3 


10. ± 2 + 2 ; 8 


10. ± 3.1 2 ; 




1.0- 


1.2 


2.4 ± 0.7±g;5 


5.± 1 + J; 4 


9.± 2 + 2 ; 





Table 5: The values of the parameter £> obtained by fitting the differential cross section to the form 
da/dp 2 oc exp(—B ■ p 2 ), together with fit errors. Because of limited statistics, no value is obtained for A 
particles produced on the W target. 





B (GcV/c)- 2 


C 


Al 


Ti 


W 




3.3 ±0.2 


3.2 ±0.3 


3.0 ±0.2 


3.3 ±0.3 


A 


2.3 ±0.3 


2.0 ±0.3 


1.7 ±0.3 


1.3 ±0.5 


A 


2.2 ±0.3 


2.0 ±0.3 


1.8 ±0.3 
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Figure 2: The Armenteros-Podolanski plot for the candidates: the transverse momentum p t of the 
oppositely charged decay products vs. their asymmetry in longitudinal momenta pf^. All momenta are 
relative to the V° line-of- flight. Background from 7 — > e+e - conversions populates the region below p t — 
0.015 GcV/c. 
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Figure 3: The invariant mass distributions for K s , A, and A particles summed over all runs with different 
targets. The size of the invariant mass bins (Am) is 3.0 MeV/c 2 for the K s distribution and 1.5 MeV/c 2 
otherwise. 



10" 



10' 



10" 



10' 









o pN =9.0+/-0.3 (mb) 




- 1 


A 




o pN =4.3+/-0.2 (mb) 






A 


o pN =1.4+/-0.1 (mb) 




10 




10 2 

A 





Figure 4: The V° total production cross sections a p A as a function of the atomic mass A of the target 
material. The solid lines show fits within the Glauber model, which yield the indicated production cross 
sections per nucleon (tr p Ar). The dashed lines are fits to the form ct p a oc A a ; the resultant a values are listed 
in Table |. 
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Figure 5: The differential cross section da p A/dp1 for K s , A, and A production on indicated target materials 
as obtained by extrapolation to the full xf range [—1, +1 ]. 
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Figure 6: Top: The ratio of A and A particle yields at mid-rapidity for pp (circles), pA (triangles), and AA 
(squares) collisions as a function of the nucleon-nucleon center-of-mass energy ^s^n- The black triangle 
denotes the average of the HERA-B results; references for the open symbols are given in the text. Bottom: 
The ratio <1(j p a(^) I 'da p A(A) determined at xf ~ —0.06 for the indicated targets; the average of the four 
values and the standard deviation are indicated by the solid and dashed lines, respectively. 




Figure 7: The total cross sections per nucleon cf p n for the production of Kg, A, and A particles as a function 
of s, the square of the center-of-mass energy. Black symbols denote the results from HERA-B, open symbols 
those from g (| |, |, [h], || || 0, [t| [It], ElI El |(| |J ||, || || and refs. therein. The curves 
are calculated using the fit functions reported in Ref. pW. 



